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The phonon density of states (DOS) and magnetic excitation spectrum of polycrystalline BiFeOa 
were measured for temperatures 200 <T< 750 K , using inelastic neutron scattering (INS). Our re- 
sults indicate that the magnetic spectrum of BiFeOs closely resembles that of similar Fe perovskites, 
such as LaFeOs, despite the cycloid modulation in BiFeOs- We do not find any evidence for a spin 
gap. A strong T-dependence of the phonon DOS was found, with a marked broadening of the 
whole spectrum, providing evidence of strong anharmonicity. This anharmonicity is corroborated 
by large- amplitude motions of Bi and O ions observed with neutron diffraction. A clear anomaly is 
seen in the T dependence of Bi-dominated modes across the Neel transition. These results highlight 
the importance of spin-phonon coupling in this material. 

PACS numbers: 63.20.kk, 75.30.Ds, 75.85.+t, 78.70.Nx 



I. INTRODUCTION 

Multiferroic materials exhibiting a strong magneto- 
electric coupling are of great interest for potential appli- 
cations in spintronic devices and actuator systems [TJ |2] . 
BiFeOa (BFO) is one of the few known systems exhibit- 
ing simultaneous magnetic and ferroelectric ordering at 
T > 300 K, and as such is a strong candidate for ap- 
plications [T]-[3]. BFO crystallizes in a rhombohedrally- 
distorted perovskite structure (space group R3c) |4HZ1- 
The Bi lone-pair is thought to be responsible for the off- 
centering of Bi atoms, which induces the ferroelectricity, 
with a high Curie temperature Tc ^ HOOK [2^. The Fe 
ions, inside oxygen octahedra, carry large magnetic mo- 
ments 4/iB [2 , and order antiferromagnetically (AF) 
below the Neel temperature, Tn — 640 K, with some 
canting of the spins, and a long-period cycloid modu- 
lation [U [HI H]. Unravelling the behavior of phonons 
and magnons, and their interactions, is crucial to un- 
derstanding and controlling multiferroic properties [l]. 
Phonons couple to the ferroelectric order, and magnons 
to the magnetic order, and it is expected that phonons 
and magnons strongly interact in a system exhibiting 
simultaneous ferroelectric and antiferromagnetic order, 
such as BFO [1^. This interaction also gives rise to mix- 
ing of the excitations, resulting in electromagnons, for 
example [D H]. To our knowledge, there are currently 
no reported experimental data of the full magnon and 
phonon spectra in BFO, however. Multiple Raman mea- 
surements have been performed, but these only probe 
modes at small wavevectors q ^ p!QH2Q] . Here, we 
report the first INS measurement of the phonon DOS 
in polycrystalline BiFeOs, as a function of T, as well 
as more detailed measurements of the magnon spectrum 
than previously reported [21 . From these data, we ex- 
tract the exchange coupling constant of BFO, and we 
identify a strong anharmonicity of the phonons, provid- 
ing evidence for strong spin-phonon coupling. 



II. NEUTRON DIFFRACTION 

A stoichiometric mixture of Bi203 (99.99%, Aldrich) 
and Fe203 (99.99%, Aldrich) was ball-milled for 10 hours 
[22] . The resulting powder was hot-pressed at 900° C for 
5 min in a half-inch diameter graphite die, with a 2 ton 
force applied, and a heating rate of 300°C/min. The 
pressed pellets were annealed in vacuum at 750 K for 24 h. 
The total sample weight was about 15 g. 

Neutron diffraction measurements were performed us- 
ing the POWGEN time-of- flight diffractometer at the 
Spallation Neutron Source (SNS), at Oak Ridge Na- 
tional Laboratory. The sample was placed inside a 
thin-wall vanadium container, and loaded in a radia- 
tive vacuum furnace. Data were collected at T = 
300, 473, 573, 623, 723, 773 K. The time-of-flight diffrac- 
tometer uses a broad incident spectrum of neutrons, and 
was configured with a center- wavelength A = 1.066 A. 
The diffraction data were refined with OS AS [23l, us- 
ing the R3c space-group. The fits indicated good crys- 
tallinity and good sample purity with ~98% BiFeOs and 
~2% of a secondary phase, indexed as Fe203, at all tem- 
peratures. For T < 640 K, the data were refined with a 
G-type antiferromagnetic order, without cycloid modula- 
tion. Since the inclusion of the G-type AF order did not 
change the refinements significantly, it is likely that the 
further incorporation of the cycloid structure would only 
have a minimal effect on our results. Results are sum- 
marized in Table [H The diffraction data and Rietveld fits 
for T = 300 K and 773 K are shown in Fig. [l] (for one of 
three detector banks). We observe a constant intensity 
ratio for (104) and (110) reflections at all temperatures 
measured, and thus our data do not support the R3c- 
R3m transition reported in Ref. [25 . 

The refined lattice parameters and atomic positions 
are in good agreement with prior reports [5 , 6 . The tem- 
perature dependences of the lattice constants and mean 
square thermal displacements (squares of quantities in 
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TABLE I: Results of Rietveld refinements of time-of- flight neutron diffraction data for BiFeOs (RSc). 
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Table [Tl) are shown in Fig. [2] The mean-square displace- 
ments were refined with an anisotropic harmonic model 
for all atoms, which assumes Gaussian probability dis- 
tributions for atom positions. The atomic mean-square 
displacements are largest for Bi, followed by O. Our re- 
sults for thermal displacements are also in good agree- 
ment with Ref. [6 , although we find displacements that 
are larger for Bi than for Fe, both along the (trigonal) c- 
axis and in the (basal) a, b plane. For Bi and Fe vibration 
modes, the behavior of {u'^) is linear in T, as expected 
in the high-T regime of an harmonic oscillator (from the 
theoretical partial phonon DOS reported in [36^, the av- 
erage phonon energies for Bi and Fe vibrations are 12 and 
28meV, respectively equivalent to 140 K and 320 K). In 
this regime, the amplitude of vibrations does not depend 
on the mass, but only on an effective force-constant, 
{v?) (xT/K [24J. Linear fits to the Bi and Fe data indi- 
cate that the effective K for Bi motions in the (a, b) plane 
is about half of the force-constant for Bi motions along c, 
which itself is comparable to those for either types of Fe 
motions. We note that our fits did not use anharmonic 
displacement models, which could slightly affect the re- 
sults. The results for {v?) of oxygen atoms in Fig. [2]-c 
shows a departure from linearity around 600 K for the 
short and long axes of the thermal ellipsoids. This effect 
may be related to the magnetic transition at Tn = 640 K, 
although there could also be effects of phonon thermal 
occupations, since the average energy of O vibrations is 
45meV, corresponding to 520 K [36 . 



III. INELASTIC NEUTRON SCATTERING 

INS spectra were measured using the ARCS direct- 
geometry time-of-flight chopper spectrometer at the SNS 
[26] . In the ARCS measurements, the sample was en- 
cased in a 12 mm diameter, thin-walled Al can, and 
mounted inside a low-background furnace for measure- 
ments at T = 300, 470, 570, 690, 750 K. All measurements 
were performed under high vacuum. An incident neutron 
energy Ei = 110 meV was used and the energy resolution 
(full width at half max.) was ~3 meV at 80 meV neutron 
energy loss, increasing to '^7meV at the elastic line. 

Additional INS measurements were performed with the 
Disk Chopper Spectrometer (DCS) at the NIST Cen- 
ter for Neutron Research, with incident energy Ei = 



3.55 meV, in up-scattering mode (excitation annihila- 
tion) [27]. In these conditions, the energy resolution was 
'^0.12 meV at the elastic line, increasing to ~1.2meV at 
25meV neutron energy gain. For DCS measurements, 
the sample was encased in the same Al can as in the 
ARCS measurements. The empty Al sample container 
was measured in identical conditions to the sample at all 
temperatures, and subtracted from the data. DCS mea- 
surements at T = 200, 300, 470, 570 K were performed in 
a high-temperature He refrigerator, and measurements 
at T = 570, 670, 720 K were performed in a radiative fur- 
nace. 



A. Magnetic Spectrum 

Figure [3] shows the orient at ion- averaged scattering 
function, S{Q^E)^ obtained with ARCS, as a function 
of temperature {Q and E are the momentum and energy 
transfer to the sample, respectively). At T = 300 K, the 
data for Q < 4A~^ clearly show steep spin- waves, em- 
anating from strong magnetic Bragg peaks at Q = 1.37 
and 2.62A~^, and extending to ^ ~ 70meV. This range 
of energies overlaps with much of the phonon spectrum 
(see below). The intensity of the magnetic Bragg peaks 
decreases with increasing T, vanishing for T > 670 K, 
in good agreement with the reported Tn = 640 K. The 
magnetic scattering nearly vanishes for Q > 6A~^, ow- 
ing to the magnetic form factor of Fe^+ ions. The high- 
E, optical part of the spin-waves is strongly damped for 
T = 470, 570 K, well below Tn. The \ow-E part of the 
spin- waves is also clearly seen in the DCS data in Fig. |4[ 
a,b, corresponding to sharp vertical streaks emanating 
from the magnetic Bragg peaks. The acoustic part of the 
spin-waves shows some Q-broadening with increasing T 
below Tn- Magnetic correlations remain for T > Tn, 
but are much broader (Figs. [3p|4]-c,d). We note that the 
orientation-averaged spin-waves show a strong similarity 
between BiFeOa and other AFeOs perovskites, such as 
ErFeOa, LaFeOa and YFeOs [MSffi. 

An important question is whether low-energy spin- 
waves are present in BFO, and whether an energy gap 
exists, potentially associated with an anisotropy in the 
exchange coupling, or single-ion anisotropy. Figure [Sj-b 
shows the magnetic scattering intensity, Sms,g{E)^ mea- 
sured with DCS, integrated over the spin- wave dispersing 
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FIG. 1: Neutron dffraction patterns (markers) and Rietveld 
refinements (red line) for BiFeOs at 300 K (a) and 773 K 
(b). Blues curves are the bottom of each panel are difference 
curves. Upper and lower tick marks are the reflection posi- 
tions for the BiFeOa and Fe2 03 phases, respectively. Peak 
labels are for the BiFeOs phase. The refinements were done 
with the R3c space-group. The data for T < 640 K were 
refined with a G-type antiferromagnetic structure (without 
cycloid modulation). 
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FIG. 2: Results of Rietveld refinements (space- group R3c) 
for lattice parameters and anisotropic atomic mean-square 
displacements, from POWGEN neutron diffraction data. 
Straight lines in (b) are fits to the data, while straight lines 
in (c) are guides to the eye. 




from Q = 1.37A~^, and compared to the phonon back- 
ground on either side. This figure clearly shows that the 
-Smagl^) intensity from the spin- wave persists down to 
|£^| 0.3 meV, where it merges with the elastic scatter- 
ing signal. Thus, we can estimate an upper-bound for any 
magnetic anisotropy gap, Eg < 0.3 meV, if any such gap 
exists. This is at odds with reports of a gap, Eg = 6 meV, 
derived from modeling the low-T specific heat Cp [31 . 

Although the powder average of excitations does not 
allow us to determine if multiple spin-wave branches ex- 
ist within the magnon DOS, we can compare our results 
with recent Raman scattering studies that reported mul- 
tiple electromagnon modes in BFO between 1.5 and 7.5 
meV [12] [13]. We do not see these excitations in the 
powder spectrum at either 200 or 300 K (see Fig. [sj-a). 
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FIG. 3: S{Q, E) for BiFeOs at different temperatures, mea- 
sured using ARCS (logarithmic intensity). 
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FIG. 4: S{Q,E) for BiFeOs at different temperatures, mea- 
sured using DCS (logarithmic intensity). 




0.0 1 ' ^ ' ^ ' ^ ' 1 

-8 -6 -4 -2 

Energy (meV) 




Energy (meV) 

FIG. 5: (a) Integrated INS intensity (DCS, 300 K) from 
spin-wave, 1.3 ^ Q ^ 1.45 A~^, compared with background 
phonon signal (1.15 ^ Q ^ 1.3 A"^ and 1.45 ^ Q ^ 1.6 A"^). 
(b) Magnetic excitation spectra obtained from INS data 
(ARCS), integrated over 2 ^ Q ^ 6 A~^, corrected for phonon 
scattering (6 ^ Q ^ 10 A~^, scaled). The thick line is the 
magnon spectrum calculation (see text). Error bars indicate 
one standard deviation. 



although high-resolution INS measurements on single- 
crystals may be needed to observe such effects. Our re- 
sults for the magnetic spectrum show a single maximum, 
around 65meV at 300 K, see Fig. [5}b. This is at odds 
with the magnetic spectrum reported in Ref. [2T1, which 
showed additional maxima around 30 and 55 meV. How- 
ever, the extra peaks in [21 are likely due to peaks in the 
phonon DOS at these energies (see below). We use the 
magnetic DOS to estimate the exchange coupling, with 
a simple spin-wave model for a collinear Heisenberg G- 
type antiferromagnet in an undistorted perovskite struc- 
ture. A similar model successfully described the spin 
waves in the orthoferrite compounds ErFeOs, TmFeOa 
[28] . as well as LaFeOa and YFeOa EH [30], which have 
close magnetic structures. This does not capture pos- 
sible effects from the spiral spin structure in BFO, but 
these are expected to be limited, owing to the long pe- 
riod of the modulation. Within this model, two exchange 
constants, J and J , describe a gapless spin- wave disper- 
sion, with acoustic and optic modes that meet at the 
magnetic zone boundaries. This behavior can be seen 
in the large patch of scattering intensity in the ARCS 
300 K data near 65 meV (Fig. pl-a). The J[J ] exchange 
constant corresponds to coupled moments with spin-spin 
distances of 3.968[5.613±0.025] A. Assuming S = 5/2 
Pg3+ ]XLoments, and comparing the maximum energy of 
the measured and model spin-wave spectra, we are able 
to place limits on the values of J and J^ We find that 
J = l.Gtoi and r = -0.25t^;?^ where there is a lin- 
ear dependence on these parameters within these bounds 
J' = 0.90(2) - J0.40(l). The calculated magnetic DOS 
for J = 1.6 and J' = —0.253 meV agrees well with the 
phonon subtracted T = 300 K data shown in Fig. |5]-b. 

There is a clear softening and broadening of the spin- 
wave spectrum with increasing temperature from 300 K 
to 570 K. Magnon- magnon and magnon-phonon interac- 
tions are likely both responsible for this softening [13 , 32j. 
It is possible that the strong softening of 6'mag(^) in this 
range is related to the an anomalous magnetization be- 
low Tn [31 . Above Tn, we observe Lorentzian scattering 
intensity centered ai E = meV, typical of paramagnetic 
behavior (the dip below 20meV is a result of imperfect 
phonon subtraction). 



B. Phonons 

In both Figs. [3] and [4] the horizontal bands of intensity 
increasing as Q cire orientation-averaged phonon disper- 
sions. In the DCS data, three main horizontal bands are 
clearly observed at|£^|~6.5,8,ll meV, corresponding to 
the top of acoustic phonon branches and low-£^ optical 
branches, which mainly involve Bi vibrations (see below). 
The acoustic branches are seen dispersing out of a nuclear 
Bragg peak at Q — 2.25A~^. The phonon cutoff corre- 
sponding to the top of oxygen-dominated optic branches 
is seen at E 85 meV in Fig. [3]-a. A strong broadening 
of the phonon modes with increasing temperature can be 
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FIG. 6: (a) Generalized phonon DOS obtained from ARCS 
data, (b) Low-^ part of generalized DOS from DCS data, 
showing strong broadening of Bi-dominated modes. Curves 
for different temperatures are vertically offset for clarity. 



seen in both Figs. |3] and [4j This broadening indicates a 
strong damping of phonons, over the full spectrum. The 
T range over which this occurs is in agreement with prior 
Raman measurements [14J, and points to a spin-phonon 
coupling effect. 

The S{Q^E) data were analyzed to extract the gen- 
eralized phonon DOS, g{E), in the incoherent scattering 
approximation. The data from ARCS were integrated 
over 6 ^ Q ^ 10 A- \ which minimizes any contribu- 
tion from magnetic scattering, and the elastic peak was 
subtracted, and extrapolated with a quadratic E depen- 
dence for £^ < 5 meV. A correction for multiphonon scat- 
tering was performed at all T [33l|34j. The gDOS from 
the DCS data was obtained by integrating over the full 
range of Q (no multiphonon correction) [35]. The mea- 
sured signal from the empty container was much weaker 
than from the sample, and was easily subtracted. The 
results are shown in Fig.jeja) for the full gDOS (ARCS) 
and panel (b) for E < 20meV (DCS). Although the DCS 
data include a magnetic contribution, this effect is lim- 
ited, and the DOSs form ARCS and DCS are in excellent 
agreement, considering the difference in instrument res- 
olution (see Fig. [t]). The coarser energy resolution in 
ARCS data washes out the three Bi-dominated peaks 
ai E ^ 15meV, but the two DOS curves are other- 
wise very similar. In BFO, the different elements have 
different ratios of cross-section over mass, cr/M, result- 
ing in a weighted phonon DOS (gDOS). The values of 
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FIG. 7: Generalized phonon DOS of BiFeOs measured 
with INS (ARCS and DCS) at 300 K, compared with first- 
principles calculation of Wang et al. |36j. 



cr/M for (Bi, Fe, O) are (0.044, 0.208, 0.265), in units of 
barns/amu, respectively. Thus, the modes involving pri- 
marily Bi motions are under-emphasized in S{Q^E) and 
g{E) (but there is relatively little weighting of Fe modes 
compared to O modes). 

While the energy-range of the spectrum is comparable 
with other Fe perovskites [30], what is striking is the se- 
vere broadening of the spectrum with increasing T. The 
gDOS measured at 300 K is in good agreement overall 
with the first-principles calculation of Wang et al. within 
spin-polarized DFT+U [36^ , as can be seen in Fig.[7| This 
agreement allows for a clear identification of the main fea- 
tures in the DOS. The three peaks at :^ 6.5, 8, 11 meV 
(Fig. |6jb)) arise from the top of acoustic branches and 
the lowest-£^ optic modes, and they are dominated by Bi 
motions. The lower two Bi peaks (6.5 meV and 8meV) 
are softer than predicted by DFT by about 10% [3^. 
These Bi modes are responsible for the peak in Cp/T^ 
around 25 K reported in Ref. [31]. The phonon spec- 
trum for E > 40meV is mainly comprised of oxygen 
vibrations, and is stiffer in the measured DOS than in 
the DFT calculation [36^. While the agreement is gen- 
erally good between the DFT calculation of Wang et al. 
and the phonon DOS measured at low temperature, the 
T dependence of the DOS is strongly affected by anhar- 
monicity and spin-phonon coupling, as we discuss in the 
next section. 



IV. DISCUSSION 

As already pointed out above, the phonon DOS ex- 
hibits a severe broadening with increasing T, which af- 
fects the whole spectrum. In particular, the broaden- 
ing of Bi modes at low-£^ is obvious in the DCS data, 
shown in Fig. [6[b. It is particularly strong between 
470 K and 570 K, with little additional broadening ob- 
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FIG. 8: (a) Centroids of the three low-energy peaks in DOS 
around 6.5, 8.5, llmeV (resp. Ei, E2, as a function 

of temperature, (b) Relative change in these energies with 
respect to their value at 200 K. The vertical dashed line at 
T = 640 K denotes the Neel temperature. 



served further above Tn- This indicates a coupling be- 
tween anharmonicity and the loss of the AF order. These 
results are consistent with the reported observations of 
strong broadening of Raman phonon modes [14 . The 
oxygen modes, dominating the DOS between 40 meV and 
85 meV, are also strongly broadened in this T-range. 

We now analyze in more detail the temperature depen- 
dence of the three peaks at 6.5, 8.5, llmeV (resp. 
£'2, Es). The peaks were fitted with Gaussians, and the 
resulting peak centers are plotted as a function of tem- 
perature in Fig. [8] As may be seen on this figure, all 
three peaks undergo a pronounced but gradual softening 
between 200 K and 570 K, besides the strong broadening. 
However, this softening stops around T^. This is com- 
patible with the recently reported behavior of Ai Raman 
modes [14]. Es actually appears to stiffen above Tn, but 
this may partly be due to the influence of broadening 
and softening of modes dit E ^ 13 meV, causing a skewed 
contribution to the Es peak. 

The change in phonon softening behavior around Tn 
is further evidence of the influence of spin-phonon cou- 
pling. The softening of phonon modes with increasing 
T, as observed here for T ^ 570 K, can generally be re- 
lated to the thermal expansion of the system through the 
quasiharmonic approximation d\nE = —jd\nV^ with 7 
the Griineisen parameter [37 . We determined the rela- 



tive change in volume to be about 1.4% between 200 K 
and 570 K from our diffraction measurements (this was 
linearly extrapolated over the range 200 — 570 K, since 
our diffraction data were limited to T ^ 300 K). The av- 
erage relative decrease of Ei, £^2, and £3 is — 4.2 ± 0.7% 
over the same T range, yielding a Griineisen parameter 
7 = 3.0 ± 0.5. Such a large value of 7 is a further cor- 
roboration of the anharmonicity of these Bi-dominated 
modes. The lack of softening above 570 K indicates that 
an increase in interatomic force-constants associated with 
the loss of magnetic order compensates for the effect of 
thermal expansion. 

We note that Bi and O atoms undergo large ampli- 
tude vibrations, according to both our measurements 
and reports of others [5l |6] . Since the Bi and O modes 
are sharp at T ^ 300 K, the large displacements ob- 
served in diffraction data are dynamic in nature, rather 
than associated with static disorder. These large ampli- 
tudes of vibration for Bi and O are related to the anhar- 
monic scattering of phonons, which leads to the broad- 
ening and softening of features in the DOS. We have also 
performed measurements of the Fe-partial phonon DOS 
with nuclear-resonant inelastic x-ray scattering (NRIXS) 
on ^^Fe-enriched samples, and observed a more limited 
broadening of Fe modes, in agreement with the smaller 
thermal displacements of these atoms [38 . 

We suggest that the large thermal displacements and 
anharmonicity of Bi and O modes lead to structural 
fluctuations, such as variations in Fe-O-Fe bond lengths 
and bond angles (tuning the superexchange interaction) 
through tilts and rotations of FeOe octahedra, that could 
lead to fluctuations in magnetic coupling. The magni- 
tude of O thermal motions perpendicular to Fe-0 bonds 
actually leads to fluctuations in the Fe-O-Fe bond angle 
that are larger (c^ 6 — 10°) than the variation of average 
angle with T. Reciprocally, the loss of magnetic order 
induces a change in interatomic force-constants, stiffen- 
ing the Bi vibration modes at low E. The motion of Bi 
atoms is also directly related to the ferroelectricity. The 
large-amplitude structural fluctuations could thus lead 
to steric effects between Bi motions and the rotations of 
oxygen octahedra, yielding a complex coupling between 
ferroelectric and AF magnetic orders. 



V. SUMMARY 

We have systematically investigated the temperature 
dependence of the magnetic excitation spectrum and 
phonon density of states of BiFeOa over the range 200 ^ 
T ^ 750i^, using inelastic neutron scattering. In ad- 
dition, we performed neutron diffraction measurements 
and refined the lattice parameters and thermal displace- 
ment parameters over 300 ^ T ^ 770K. We separated 
the magnon and phonon contributions in the S{Q^E) 
data, and observed a strong resemblance of the magnon 
spectrum with that of related compounds LaFeOa and 
YFeOa. The magnon spectrum was fit satisfactorily with 
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a simple collinear Heisenberg model for a G-type antifer- 
romagnet, indicating the limited role of the cycloid on 
the spin dynamics, as expected from the long period of 
the cycloid modulation. The phonon DOS obtained from 
the high-Q part of the S{Q^ E) data is in good agreement 
at low temperatures with the first-principles calculation 
of Wang et al. [36 . However, the phonon DOS shows 
a strong temperature dependence, with in particular a 
pronounced broadening. Also, both the softening and 
broadening of features in the DOS correlate with the loss 
of antiferromagnetic order around T/v = 640 K, indicat- 
ing the presence of significant spin-phonon coupling, in 
agreement with recently reported Raman measurements 
[M]. The potential infiuence of large atomic displace- 
ments on the modulation of the superexchange interac- 
tion, and the concomitant effect of the change in force- 



constants from the loss of magnetic order were pointed 
out. 
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